Emission angle distribution and flavor transformation 

of supernova neutrinos 



Wei Liao 

Institute of Modern Physics, East China University of Science and Technology, 
P.O. Box 532, 130 Meilong Road, Shanghai 200237, P.R. China 

Center for High Energy Physics, Peking University, Beijing 100871, P. R. China 

Abstract 

Using moment equations we analyze collective flavor transformation of supernova 
neutrinos. We study the convergence of moment equations and find that numerical 
results using a few moment converge quite fast. We study effects of emission angle 
distribution of neutrinos on neutrino sphere. We study scaling law of the amplitude 
of neutrino self-interaction Hamiltonian and find that it depends on model of emis- 
sion angle distribution of neutrinos. Dependence of neutrino oscillation on different 
models of emission angle distribution is studied. 

PACS: 14.60.Pq, 97.60Bw 

1 Introduction 

Flavor transformation of neutrinos in core-collapse supernova is one of the important 
remaining problems in neutrino physics. This problem has been investigated by many 
researchers [D El El l El El El S El EDI El [iSl El [11 [151 [161 ^ It 
is realized that neutrino density above neutrino sphere in supernova can be so large that 
neutrino-neutrino refraction can dominate flavor transformation of neutrinos. Research 
on the effect of neutrino-neutrino refraction is difficult because it is caused by neutrino 
self-interaction and is of non-linear nature. Complete numerical analysis use discrete set 
of energy bins and angle bins of neutrinos. Evolutions of million equations have to be 
followed. It is very complicated. 

In a recent work we derived a set of moment equations describing the transport and 
ffavor transformation of neutrinos in supernova |23]. Distribution of neutrinos over angle 
9p, the angle of neutrino direction intersecting with radial direction in supernova, is en- 
coded in moments of density matrix. The equation of neutrino is expanded using these 
moments. Instead of using a large number of angle bins we just need a few moments to do 
numerical study. It is shown that numerical works can be simplified by about two orders 
of magnitude in comparison with multi-angle simulation. Moreover, this formulation of 
neutrino in supernova also offers us a way to study the effect of emission angle distribution 
of neutrinos on the transport and ffavor transformation of neutrinos. 
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In this article we analyze the effect of emission angle distribution of neutrinos. In 
section [2] we make a quick review on the moment equations. In section [3] we analyze the 
scaling behavior of the strengths of moments in different models of neutrino emission. 
We check convergence property of moment equations in the analysis of the strengths 
of moments. In section H] we analyze effect of different models of neutrino emission on 
collective neutrino oscillation. We summarize in section [51 



2 Moment equations 

In Ref. p3j we introduced moments of pp{t, r), density matrix for neutrinos at given time 
t and radius r: 

Pkit,r,\p\) = j da^ {l-cos9p)' pj,{t,r), k = 0,1,2- ■■ , (1) 

where dp is the angle of neutrino direction intersecting with the radial direction, as shown 
in Fig [H Similarly we introduced pk for anti-neutrinos. We also introduced re-scaled 
moments 

Pi = z'^'^'^p,, (2) 

where 

z = r/ro, (3) 

tq is the radius of neutrino sphere. Similarly we introduced p'^ for anti-neutrinos. In Fig. 
[Done can see clearly 



sin 6p = — sin 6po (4) 



It is easy to see 



1 - cos Op = (^rl/r'^ j + ^l-^sin^^^o ) • (5) 

It scales approximately as r^^. Together with the scaling behavior of the zeroth moment 
a geometric factor ^-2(fc+i) fQ^j^j;^ for moment pk- The factor ^^{fc+i) jg introduced in 
Eq. ([2]) to compensate this geometric scaling factor. 

Using some approximations we arrive at the following set of moment equations 

^ =-r,^Ql-i[HA,p'ul k = 0,l,---,N (6) 
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Figure 1: Geometric picture of angles of tlie neutrino momentum intersect with f. 



where > 1 is an integer 



N 



Ql = z 



.2k 



(7) 



l=k+l 



Ha = Ho + V2Gf{L + z-^Di), 




I 



dE 



[p;(r,E)-p;(r,E)]. 



(9) 



(27r)3 



Eq. ([6]) is a set of truncated moment equations in P/v approximation for which p'^ = 
(Pfc ~ 0) been set for k > N. Q\f = 0. Hq is the Hamihonian for vacuum oscillation, 
L = diag{ne, n^^rir} in the flavor base is the matter term given by charged lepton number 
densities Ue^n^r- Gp is the Fermi constant. Equation for pk is similar except replacing Hq 



A few points concerning moment equations are as follows: a) Physical observables are 
described by po and pi. Integration of E'^Tr[po\ over energy gives the neutrino density 
and integration of E'^Tr[pQ — pi] gives the neutrino flux; b) Emission angle distribution of 
neutrinos on neutrino sphere is described by moments pk and their effect in the neutrino 
flavor transformation can be systematically studied; c) The strength of p^, Tr[py, is 
modified by Ql term and does not change if this term is neglected; d) The scaling law of 
the self-interaction Hamiltonian is no longer when N > 1 and is modified by higher 
moments. Precise scaling behavior should depend on the model of neutrino emission. 



by -Ho. 
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Figure 2: (color online) Fo and Fi in Model I. 

3 Emission angle distribution and scaling law of mo- 
ments 

In this section we analyze the strengths of zeroth and first moments, that is Tr[po,i], in 
different models of neutrino emission. This analysis can tell us a lot on how strong neu- 
trino self-interaction is. It can also tell us a lot on the convergence property of moment 
equations. This is because po and pi are the most important quantities in our problem. 
Physical observables are given by po and pi. When neutrino self-interaction gives domi- 
nant contribution flavor transformation of neutrinos is controlled by Di which is directly 
related to pi. Effects of higher moments on Tr[po,i] tell us how large higher moments 
affect the flavor transformation of neutrinos. 

We consider three models of neutrino emission on neutrino sphere. 

Model I, neutrino is uniformly emitted with respect to the emission angle 6pQ and 

Pkit,ro) = j^^poit,ro) (10) 
Model II, emission angle distribution of neutrinos is proportional to cos^po and 

P.(t,ro) = ^^^^^^Po(t,ro), (11) 
Model III, emission angle distribution of neutrinos is proportional to (1— cos 6*^0) cos6'po 
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Figure 3: (color online) Fq and Fi in Model II. 



and 

6 

P.(t,ro) = ^^^^y^^Po(t,ro). (12) 
The evolution of Tr[pk\ is simple and is obtained by taking the trace of Eq. ([6]): 

^ = -o '0. (13) 

where 

N 

Qu = z'^ ^(/ + l).-(^'+i)Tr[p;]. (14) 

l=k+l 

Qn = 0. The second term in ([6]) does not contribute to Tr[py. 
We do numerical analysis for Fq and Fi. 

Fo = Tr[p',ir)]/Tr[p',iro)], F, = Tr[p[ir)]/Tr[p',iro)] (15) 

Fq^i are Tr[pQ^] relative to Tr[p'(^] at r = tq. In our numerical analysis we work in two 
flavor system of (i/g, u^)- We choose L^^ = Lp^ = L^^ = Lp^ = 3. x 10^^ erg/s. The initial 
energy spectrum of neutrino is given by the Fermi-Dirac distribution 
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Figure 4: (color online) Fq and Fi in Model III. 



where x = E/Ty and N2 is the normahzation factor. Parameters of four types of neutrinos 
and anti-neutrinos are chosen as: Ty^ = 2.76 MeV, Tp^ = 4.01 MeV, T^^ = Tp^ = 6.26 
MeV. fiy^ = fip^ = /i^^ = /ip^ = 3. 

In Fig. [2], [3] and H] we show resuhs in Model I, II and III separetely. A number of 
characteristics can be read out in these figures: 

i) In Pi approximation Fi keeps as a constant. This is because Qi is set to zero in 
this approximation. 

ii) In P2 approximation Fi is modified. Results of Fq do not agree with those in Pi 
approximation. 

iii) In P3 approximation results of Fq become close to those in P2 approximations. 
Results of Pi do not agree with those in P2 approximation. This is because in P2 approxi- 
mation Pi becomes corrected by Tr[p'2\ but Tr[p2] is still a constant. In P3 approximation 
Tr[p2] is also corrected and its contribution to Pi is modified. 

iv) Results of Model II and III in P4 and Pg approximations agree perfectly for both 
Fq and Pi. 

v) Results of Model I in P4 and Pq approximations are in good agreement for Pq. For 
Pi there are still some small differences. 

A few comments are as follows: 

a) Results in Model II and Model III converge faster than the results in Model I. This 
is in agreement with the observation that higher moments in Model II and Model III are 
more suppressed than those in Model I. Hence Model II and Model III should have better 
convergence properties. 
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b) Value of Fq at large radius can be understood using flux conservation. The flux of 
neutrino is given by 

Tr[po - pi] = z-'Tr[p', - z-'p[] (17) 

The flux, as it should be, scales as z~'^{ or r~^) in stationary approximation. So Fq — z~'^Fi 
is a conserved quantity. At large r this quantity approaches to Fq. On the other hand its 
initial value can be read out directly from the models of neutrino emission. Using Eqs. 
(Ho]), dH]) and (ini) we find that at large r 

Fo ^ ^, in Model I (18) 
2 

Fo ^ -, in Model II (19) 
3 

Fo ^ ^, in Model III (20) 

These values are in agreement with the plots in Figs. [21 [3] and HI 

c) The scaling behavior of Fi tells us that in Pjv approximation with > 1 the 
self- interaction Hamiltonian scales down faster than rg/r*^. 

d) Numerical study shows that Tr[p'i^] with k > 1 also drops down by 10^"*^ — 10^^ at 
large r. It is a further support to the point that moment equations converge quite fast. 



4 Flavor transformation 

In this section we do some analysis on flavor transformation of supernova neutrinos. We 
study the case of inverted mass hierarchy and for simplicity we neglect matter effect in 
the analysis. 

In Fig. [5] we give plots of z/g fraction versus radius. These plot are obtained by solving 
Eq. ([6]) numerically. For a small step we get 

p',{r + Ar) = -—QMr)] + e^^^ V.(r)e^^-^'- (21) 

In these plots one can see synchronized oscillation for which neutrinos of all energy 
point to the same direction in flavor space. Beyond the region of synchronized oscillation 
neutrino flavor vectors spin down which leads to neutrino flavor conversion. 

We compare numerical results of P4 approximation and of Pq approximation in models 
I and III. We flnd nice agreements between these two approximations. This shows that 
Pn approximation converge quite fast. This is in agreement with the discussion in the 
last section that the scaling law of the strength of the Hamiltonian converge quite fast. 
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Figure 5: (color online) Fraction of v^, ny^ jijiy^ + rij^^), versus radius r in in different models. 
Left in Model I; Right in Model IL | Am|J = 3 x lO^^eV^ sin^ 26^13 = 0.01. 

In Fig. [5] one can see that result of Pi approximation is quite different from that of P4 
and Pe approximations. This is also consistent with discussion in the last section. Since 
the scaling law of the Hamiltonian in Pi approximation is quite different from that in P4 g 
approximation we would expect to find difference in oscillation pattern. 

In Fig. Owe compare numerical results in models I, II and III. One can see that there 
are some differences in the oscillation pattern. This is consistent with the analysis on the 
strength of self-interaction Hamiltonian in these models. Numerical results show that at 
large r 

Pi ^ 0.126, in Model I (22) 

Pi ^ 0.134, in Model II (23) 

Pi ^ 0.151, in Model III (24) 

Since differences in Hamiltonian are not large at large radius and the differences in the 
oscillation pattern should not be large either. 

We note that the scaling law of Pi, hence the amplitude of neutrino self-interaction 
Hamiltonian, is model dependent. This dependence on model is nicely described by the 
corrections given by higher moments in moment equations. Previous researches use fixed 
scaling function for the self-interaction Hamiltonian and do not take into account the 
dependence of the scaling law on the emission angle distribution of neutrinos. As a 
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Figure 6: (color online) Fraction of Vf,, Uy^jiny^ + it-u^): versus radius in different models. 
Neutrino parameters are the same as in Fig. [5j 

comparison one can check the result using a fixed scahng function. For example, one 
can use pi = 5Vi(^o) wher e S{r) = ^7(1 + Hence Fi = 5/(1 + Vl^^ 

in model = 0.33/(1 + Vl - z^f in model II and Fi = 0.5/(1 + Vl - z^f in model 
III. At large radius this model independent scaling function gives Fi 0.125 in model I, 
Fi 0.083 in model II and Fi 0.125 in model III. Only for model I this fixed scaling 
function gives a correct result at large radius. In model II this fixed scaling function gives 
a value quite different from the value obtained using moment equations. 

5 Conclusion 

In summary we have analyzed some properties of moment equations and the flavor trans- 
formation of supernova neutrinos. We have analyzed the scaling behavior of neutrino 
density and the amplitude of self-interaction Hamiltonian of neutrinos. They are related 
to quantities po and pi. 

We analyzed the convergence of P/v approximation of moment equations. Numerical 
results show that the scaling behavior of Tr[po,i] converge for N < 10. We show that 
results of neutrino oscillation also converge quite fast. These analysis are consistent. Since 
the integration of E^pi give the self- interaction Hamiltonian the analysis on Tr[pi] tell us 
how fast the amplitude of neutrino self-interaction converge. 

We analyze neutrino flavor transformation. We find synchronized oscillation and bipo- 
lar oscillation in the oscillation pattern of supernova neutrinos. We find that oscillation 
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pattern of neutrinos converge quite fast for < 10. The Pi approximation can be used 
to make qualitative analysis but can not be used to do precise numerical study. 

We study three models of emission angle distribution of neutrinos on the oscillation 
pattern and analyze model dependence of neutrino flavor transformation on the emis- 
sion angle distribution of neutrinos. Different models of emission angle distribution can 
give different results in the scaling behavior of self-interaction Hamiltonian and in oscil- 
lation pattern of neutrinos. This model dependence is carefully taken into account in the 
correction given by higher moments in moment equations. 

Previous works on oscillation of supernova neutrinos use fixed scaling function for the 
self-interaction Hamiltonian and do not take into the fact that the scaling law can be 
different in different models of neutrino emission. Analysis on the model dependent effect 
of emission angle distribution in neutrino oscillation is not presented in previous works. 

Acknowledgement: I wish to thank Y. Z. Qian, G. RafFelt, A. Yu. Smirnov and H. 
Duan for discussions on neutrino flavor conversion in supernova. 
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